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ABSTRACT 

Fluctuations in brightness due to water vapor in the Earth's atmosphere are an 
important source of noise for ground-based instruments attempting to measure the 
anisotropy of the Cosmic Microwave Background. This paper presents a Kolmogorov 
model of atmospheric fluctuations, and derives simple expressions to predict the impact 
of fluctuations on experimental measurements for three instrument configurations: 
chopped beam, swept beam and interferometer. Data from the South Pole and from 
the Atacama Desert in Chile, two of the driest places on Earth, are used to characterize 
the fluctuations at each site. Using an interferometric instrument as an example, the 
data suggest that the South Pole is the superior site for observations of the Cosmic 
Microwave Background at degree angular scales. 

Subject headings: atmospheric effects — instrumentation: interferometers — methods: 
analytical — cosmic microwave background 

1. Introduction 

Many groups are currently engaged in measuring the level of anisotropy present in the 
cosmic microwave background (CMB) from arcminute to degree angular scales. The brightness 
temperature of the fluctuations is of order 10~^ K, requiring carefully designed experiments. 
Contamination from point sources and galactic dust emission is minimized by choosing an 
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observing frequency of between 20 and 300 GHz (1.5 cm to 1 mm), depending on the angular 

scale being investigated (Tegmark & Efstathiou 1996). Differential measurements are employed to 
minimize systematic errors: chopped beam instruments measure the difference in emission between 
two or more directions on the sky, swept beam instruments sweep a beam rapidly backwards 
and forwards, and interferometers correlate signals from two or more antennas. Some examples 
of chopped beam experiments are Python I-IV (Dragovan et al. 1994, Ruhl et al. 1995, Piatt 
et al. 1997, Kovac et al. 2000), OVRO RING5M (Leitch et al. 1998), and Tenerife (Davies et 
al. 1996). Current or recent swept beam experiments include Python V (Coble et al. 1999), 
Saskatoon (Netterfield et al. 1997), and the Mobile Anisotropy Experiment (Torbet et al. 1999). 
Interferometers currently being built include the Degree Angular Scale Interferometer (DASI) 
(Halverson et al. 1998), the Cosmic Background Interferometer (CBI) (A.C.S Readhead & S. 
Padin, personal communication), and the Very Small Array (VSA) (Jones & Scott 1998). In 
addition, the Cosmic Anisotropy Telescope (CAT) has already been operational for several years 
(Scott et al. 1996). These employ wide bandwidth, low-noise receivers to maximize sensitivity, 
and should be deployed at sites where the Earth's atmosphere does not significantly compromise 
performance. 

The atmosphere is also a source of brightness temperature variations, originating primarily 
from water molecules. Of the the three states that may be present - vapor, liquid and ice - it 
is water vapor that is most important. Most of it is contained in the troposphere with a scale 
height of ~ 2 km, and, because the water vapor is close to its condensation point, it is poorly 
mixed with the 'dry' component of the atmosphere (mostly nitrogen and oxygen). This, in 
combination with turbulence, leads to a clumpy, non-uniform distribution of water vapor in the 
troposphere. Since the water molecule has a strong dipole moment, rotational transitions couple 
strongly to millimeter-wave radiation, and water vapor is the dominant source of atmospheric 
emission (and therefore opacity) at most millimeter wavelengths. Liquid water, in the form of 
clouds, is also a source of non-uniform emission, but radiates much less per molecule. Ice is the 
least efficient radiator, since the molecules are unable to rotate. Fluctuations in temperature 
caused by turbulent mixing are also a source of brightness temperature variations, although they 
are generally nuich less significant than the water vapor contribution at millimeter wavelengths. 

It should also be noted that the high refractive index of water vapor causes an excess 
propagation delay, and a non-uniform distribution of water vapor distorts an incoming wavefront. 
This sets a 'seeing' limit on interferometric observations at millimeter wavelengths, and there is an 
ongoing effort to correct for this effect. It is this problem, which limits performance at arcsecond 
spatial resolution and high frequencies, that has driven much of the recent research into the 
distribution of water vapor (e.g. Armstrong k, Sramek 1982; Treuhaft & Lanyi 1987, Wright 1996, 
Lay 1997). 

The wavefront distortions are not significant for the low resolution experiments considered 
here. The atmospheric emission fluctuations, however, can only be distinguished from fluctuations 
in the CMB by the wind-induced motion of the atmosphere with respect to the background. 
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This paper investigates how well the two can be separated. The next section describes a model 
of atmospheric fluctuations and the responses of the different types of instrument. Section ^ 
describes how data from the Python V experiment have been used to characterize the fluctuations 
at the South Pole, and Section Q estimates the level of emission fluctuations for the Atacama 
Desert in Chile using rms path fluctuation data. Section ^ shows how the theory of Section |2| can 
be combined with the fluctuation data from each site to predict the residual noise level due to the 
atmosphere for a given instrument conflguration. 



2. Models 

In this section, we develop a model that describes the atmospheric fluctuations and their 
interaction with the instruments commonly used to measure CMB fluctuations. 

Church (1995) presented a rigorous mathematical analysis based on the autocorrelation 
function of the fluctuations. We adopt a more pictorial approach based on the power spectrum of 
the fluctuations, using an atmospheric model that differs in a fundamental way from that assumed 
by Church. 



2.1. Sky brightness and instrument response 

A given pointing on the sky can be represented by a point on the surface of a sphere. Consider 
a patch of sky with angular extent /S.9x <C 1 radian, /S.6y <C 1 radian, such that the curved surface 
is well represented by a pair of ortholinear angular coordinates {9x,0y). In addition, we will 
consider only the fine structure (angular scales <^ 1 radian) in the sky brightness distribution. 
These approximations greatly simplify the analysis; simple Fourier Transforms can be used, and 
the gradient in brightness temperature due to airmass can be ignored, without significant impact 
on the degree scales of interest. 

The Fourier Transform of the (fine structure) sky brightness distribution Ts]^y(9x,0y) is given 

by 

T,^y{ax,ay) = ^=== / / T,^yi0x,9y)e-'^^^'^^'^+"y'yUexdey. (1) 

y^AOx^fy J-Ae^/2 J-Aey/2 

The angular wavenumbers {ax,ay) have units of cycles per radian. In general Tsky(aa',ay) is a 
complex variable. Its amplitude is denoted by Tsii,y{ax,ay). The normalization in equation (|l]) is 
chosen such that the variance of the sky brightness fluctuations is given by 

I r+Ae^/2 r+A0y/2 

-iZ II rjiZ 



^sky.rms — Afl A /3 / / Tky{6x, 9y)d9xd9y (2) 

^' ^9xA9y J-Ae^/2 J-Aey/2 ^ 

f^i,y{ax,Oiy)daxday, (3) 
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where Tj^^ represents the Power Spectral Density (PSD) of the (fine structure) sky brightness 
distribution (units: temperature^ radian^). 

A radiometer has a beam pattern represented by gain function G{9'^,9y), where 9'^ and 9'y 
are angular offsets from the optical axis of the instrument. For a beam pattern localized within 
9'^ <^ l,6y <^ 1, the response to fluctuations of a given angular wavenumber is 

/+00 r+oo 
/ G{e'„ 9'y)e~^-^^"^'^+»y'yU9'j9'y. (4) 

The gain G{9'^,6y) is normalized such that the maximum value of G{ax,ay) is unity. The hat 
symbol is used instead of the tilde to distinguish between the different normalizations applied 
in equation (|^) and equation The response of the instrument to the sky is given by the 
convolution of the beam and the sky brightness. In the angular wavenumber domain, 

Tinst = TskyG. (5) 

The instrument acts as a spatial filter that is only sensitive to a particular range of angular 
wavenumbers on the sky. For example, a wide-angle beam smears out the small-scale structure, 
and the corresponding G has a narrow distribution. An ideal pencil beam responds equally to all 
angular wavenumbers with G = 1. 

Atmospheric features in the sky brightness distribution are blown across the angular 
coordinate frame, causing the output of the radiometer to vary with time. Features on the sky 
with large angular wavenumber (many cycles per radian) produce rapid variations with time 
compared to those with small angular wavenumber. Time averaging of the instrument output 
over a period tav suppresses the rapid fluctuations (large wavenumber). This relationship can be 
expressed by including a temporal filter function Q{ax,ay), i.e. 

/+00 r+oo 
/ fl-yG^Q^daxday. (6) 
-OO J —CO 

Equation (|6|) can be used to calculate the rms noise due to atmospheric emission fluctuations 
for a radiometric system with time averaging. Section 2.2 describes an approximate model for T'gky 
Section 2.3 derives the spatial filter function G for chopped beam, swept beam and interferometric 
experiments. Section 2.4 characterizes the instrument temporal response Q. Section 2.5 applies 
the results of 2.2, 2.3 and 2.4 to derive the residual level of atmospheric fluctuations for a 
measurement, using equation (^). Section 2.6 highlights the atmospheric parameter most relevant 
for degree-scale imaging of the CMB. 

2.2. Model of Atmospheric Emission Fluctuations 

We adopt the Kolmogorov model of turbulence (Tatarskii 1961), summarized briefly below. 
Turbulent energy is injected into the atmosphere on large scales from processes such as convection 
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and wind shear, and then cascades down through a series of eddies to smaller scales, until it is 
dissipated by viscous forces on size scales of order 1 mm. If energy is conserved in the cascade then 
simple dimensionality arguments can be used to show that the power spectrum of the fluctuations 
in a large 3-dimensional volume is proportional to q~^^^^, where q is the spatial wavenumber 
(units: length"^). This holds from the outer scale size Lq on which the energy is injected to the 
inner scale size Li on which it is dissipated. Tatarskii showed that this same power law applies to 
quantities that are passively entrained in the flow of air, such as the mass fraction of water vapor. 

Figure |l| shows the geometry used for this analysis. The water vapor fluctuations are present 
in a layer of thickness Ah at average altitude /lav The x-, y- and z-axes form an orthogonal set 
with the z-axis parallel to the line of sight of the observations at elevation e. 

Consider first the case of observing in the zenith direction, e = 90°. In the optically thin limit, 
the brightness temperature contribution at frequency u from a volume element of atmosphere with 
thickness dz along the line of sight and water vapor density P112O is given by 

dTs]^y{x,y,z,u) = pu2o{x,y, z)Kii2o{i^)Tphys{x,y, z)dz, (7) 

where kh2o('^) is the mass opacity function for water vapor and Tphys is the physical temperature 
of the volume element. The brightness temperature of the atmosphere looking vertically up from 
the ground is 

roc 

Tsky{x,y,i^) = Kii^oi^) Tpi,ys{z)pH2odz. (8) 





The dependence on v is considered to be implicit from now on. The Fourier transform of this 
distribution is 



Tsky{qx,qy) = ^lim 





r+y/2 


r+X/2 


VXY J 


'-y/2 J 


'-X/2 



Tsky{x, y)e-'^^'^''-''+''yyUxdy 



(9) 



where Qx and Qy are the spatial wavenumbers for the x and y directions. 



For Kolmogorov turbulence it can be shown (e.g. Lay 1997) that the Power Spectral Density 
(PSD) for the fluctuations as seen in projection is given by 

The coefficient A, and the related value A', are a measure of the turbulent intensity. (T!,ky) has 
units of m^ and the angular braces denote an average over many realizations of the random 
atmosphere. The first case in equation (|lO|), applies to spatial wavelengths smaller than 2 Ah, 
where the turbulence is considered to be isotropic in three dimensions. In the second case the 
turbulence is isotropic in the horizontal plane, but is constrained vertically to lie in a layer of 
thickness Ah. The increase of the PSD as {qx,Qy) — > (0)0) is less rapid than for the 3D case. 
Beyond the outer scale Lq the PSD should become constant. 
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The derivation of the —11/3 power law by Kolmogorov apphes in the isotropic three- 
dimensional case, well within the inner and outer scales, where the turbulence can be considered 
scale-free. There is extensive experimental evidence to support this. In previous analyses 
(e.g. Church 1995, Andreani et al. 1990), it was assumed that there were no correlations present 
in the turbulent layer on scales greater than the thickness Ah, i.e. Ah corresponded to the outer 
scale size (outer scale sizes of between 1 m and 100 m were used in Church's calculations). Data 
from atmospheric phase monitors operating at 12 GHz (Masson 1994; Holdaway 1995, Lay 1997) 
and radio telescope arrays such as the Very Large Array (Armstrong & Sramek 1982; Carilli & 
Holdaway 1997) show no evidence for such a low value of Lq; indeed, correlation is observed over 
separations in excess of 10 km. While there may be conditions in which the outer scale length 
is greatly reduced (Coulman &: Vernin 1991), the data are generally well described by a —11/3 
power law on small scales, and a —8/3 power law on large scales, with the transition occurring 
for sizes comparable to the layer thickness as described in equation (p!o|). The —8/3 power law 
for the two-dimensional regime can be derived from similar scaling arguments to the 3D case, 
but cannot be justified well on theoretical grounds, since turbulence is not strictly possible in 
a two-dimensional medium - the layer must have some vertical extent - and the model should 
therefore be regarded as somewhat empirical. This distribution of fluctuation power is shown 
schematically in Fig. |2|. 

Linear coordinates {x,y) for fluctuations in a layer at average height above the ground, /lav, 
are converted to angular coordinates in radians: 6x — x/ha.^,6y ~ y/hs.Y (still assuming e ~ 90°). 
The approximation holds when Ox ^ 1 and Oy <^ 1. Linear wavenumbers {qx,Qy) are converted to 
angular wavenumbers: Ox = Qxhsw, cty = <lyha.v Therefore 

,f2 ( Uhll\al + aiy/^ /iav/(2A/i)«(a2+a2)V2«„. 

\ skyl«x,«,j; \A'hli\al + alr^/<' « («^ + a^)^/' « /iav/(2A/i). ^ ' 

The transition between the three- and two-dimensional regimes is at approximately 
/iav/(2A/i); Oo and Oi correspond to the outer and inner scales of the turbulence, 
respectively. A factor of /i^v is included for the correct normalization of the PSD, such 
that {f^^^y{ax,ay))daxday = {f^]^y{qx,qy))dqxdqy. 

Up until this point it has been assumed that the observations are in the zenith direction. 
In the 3D regime, two factors are needed to scale to arbitrary elevation. First, the atmospheric 
fluctuation power is proportional to the path length through the layer, which scales as 1/ sin e, 
i.e. A in equation ( |ll| ) must be replaced by A/ sine. In addition, the distance to the layer of 
fluctuations has increased from /lav to /lav/ sine. Applying these substitutions to equation (pT|) 

4) 



for the 3D regime only (/iav/(2A/i sin e) <C {al + a^)^/^ <C aj) 



{TLiax,ay)) = (^)(^f\al + al)-''r^ (12) 



Ahll,\,ine)-^l^a-l^'\ (13) 
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Fig. 1. — Geometry adopted for the analysis. The generic instrument is observing at elevation e 
through a layer of fluctuations with thickness Ah at average altitude /lav that is being blown by 
wind vector w. The x-, y- and z-axis form an orthogonal set with the z-axis parallel to the line of 
sight, and the a;-axis parallel to the projected wind direction. 




Fig. 2. — Schematic representation of the distribution of atmospheric fluctuation power in the 
angular wavenumber domain. The combination ax = 10, ay = would correspond to a fluctuation 
with a periodicity of 0.1 radian (5.7°) in the .x-direction. In practice the distribution is distorted 
within (a^ + ay)^^'^ ^ 1, since the angular approximation breaks down. 
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This prescription for atmosphere fluctuations, valid for small angular scales, is used in §2.5. On 
larger scales, where the 2D regime becomes more important, the situation has a much more 
complicated dependence on layer thickness and elevation, and a numerical calculation should be 
performed (e.g. Treuhaft and Lanyi 1987; Lay 1997). 



2.3. Spatial Filtering 

This section is concerned with the derivation of the spatial filter function G'^{ax,ay) 
(eq. 1^) for three common instrument configurations: (1) chopped beam, (2) swept beam, and 
(3) interferometer. The function G{ax,ay) is the Fourier Transform of the instrument gain 
pattern, G{9x,0y) (eq. 0). We assume that the emission fluctuations are present in the far-field 
of the instrument, i.e. /lav so that this analysis does not apply to large telescopes and short 

wavelengths. 

The result for a single aperture with a circular, Gaussian beam pattern is useful for all 3 
configurations. For a beam with Full-Width-to-Half-Maximum (FWHM) power of 9\^, it can be 
shown that the spatial filter is another circular Gaussian: 

^ / ^'gg(ag. + «g) l 

G iax,ay) = expi 2b^2^ f' ^ ' 

with FWHM of 0.626'b^ 



2.3.1. Chopped beam 

A single aperture is chopped between two positions on the sky at equal elevation separated 
by angle ^chop and the output is the difference between the two signals. The power gain pattern 
G{9x, 9y) is shown schematically in Fig. with the 9x axis defined to be along the chop direction. 
The individual beams are assumed to be circular Gaussians with FWHM of ^b- 

The corresponding spatial filter G^{ax,ay)^ depicted in Fig. |^, is derived using the Fourier 
Transform relationship in equation (§), and is given by 

G'^{ax,ay) = sin^(7r6lchopax)exp < - ^^^^ ^ > . (15) 

The normalization is such that the maximum value of G{ax, cxy) is unity (Section |2.1| ). The lower 
half of Fig. 1^ shows a cross section of the spatial filter along the ax axis. The instrument has 
zero response to brightness corrugations on the sky that have ax = n'/9~^^p, (n = 0, 1, 2, . . .). It is 
also insensitive to scales much smaller than the beam size 9^, corresponding to large values of a. 
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Fig. 3. — Schematic responses of a chopping single aperture (a), a swept beam single aperture (b) 
and an interferometer (c) to far-field emission on the sky. The square of the Fourier Transform of 
each sky response is shown in (d), (e) and (f). A profile of the spatial response along the a^^-axis is 
also shown. The arrow in the upper plots indicates the projected wind vector that gives rise to the 
temporal filter (thin dark strip) shown in the lower figures. Circles represent the half-maximum 
contour of distributions that are approximately Gaussian. The shaded region in (e) represents one 
channel in the power spectrum, with an effective width Aachan that is larger than the channel 
spacing ^^eep if tapering has been applied to the swept 'slot' in (b). See section 2.3 for a full 
description of these figures. 
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2.3.2. Swept beam 

The swept beam configuration sweeps tlie beam from a single aperture back and forth across 
a strip of sky of width ^sweep (Fig. ^). The single beam is assumed to be a circular Gaussian 
with FWHM of 9^, and the Ox axis is defined by the sweep direction. It is assumed that the data 
are sampled at the Nyquist rate to generate N = 20swecp/^b samples per sweep. An FFT of these 
data generate real channels (or alternatively N/2 complex channels for a^. > 0) in the angular 
wavenumber domain spaced by ^^eep (dashed lines). The resulting spatial filter is illustrated in 
Fig. This case is different from the single chop, since there are N independently determined 
quantities, rather than one. Each channel is sensitive to corrugations on the sky with a particular 
range of wavenumber ax- In practice, data from the swept 'slot' on the sky must to be tapered 
at the ends to avoid 'ringing' effects in the transform. This gives rise to an effective channel 
width Aa that is larger than the channel separation, as indicated by the shaded region in Fig. 
The beam size again sets an upper limit on the angular wavenumber to which the instrument is 
sensitive. 

2. 3. 3. Interferometer 

Figure ^ and f show the response of an interferometer consisting of two circular apertures 
producing Gaussian beams with FWHM of ^b- The aperture centers are separated by a baseline 
of length B and the Ox axis is defined to be in the same direction. The interferometer responds 
to a range of spatial wavenumbers centered on = ±B/X (Fig ^). The larger the aperture, 
the larger the range of wavenumbers the instrument is sensitive to. It is assumed that a complex 
correlator is used to measure both the sine and cosine components on the sky; half of the total 
fluctuation power is present in each component (for the sake of clarity. Fig. ^ shows only the cosine 
component). Note that the Gaussian profile would imply that the interferometer has a finite 
response to {ox, a^) = (0, 0). In fact there must be zero response and the Gaussian approximation 
breaks down near the origin, since there is no overlap between the apertures. 

2.4. Temporal Filtering 

This section derives the form of the temporal filter Q{ax, cty) of equation (^). 

The wind vector w advects the layer containing the water vapor fluctuations in a horizontal 
direction, so that the distribution of fluctuations projected onto the x — y plane appears to move 
at speed Wxy = {wx, Wy), the component of w parallel to the x — y plane. A fluctuation component 
characterized by wavenumbers (g^;, (?j^) sampled along a line of sight parallel to the z-axis produces 
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a signal that varies in time with frequency v = WxQx + "^yQy Converting to angular wavenumber, 

sin e 

= -7 [Wxax + WyOy). (16) 

"-av 

Time-averaging the output of an instrument is equivalent to a low pass filter which rejects 
signals that are varying rapidly. Since time-averaging is equivalent to convolution of the output 
time series with a boxcar function, the equivalent frequency response is given by sinc(7rz^tav)- 
Substituting for the frequency u we obtain the temporal filter function Q: 

Q(ax,ay) = sine ^ '^^av sm e ^^^^^ ^ WyOy)^ . (17) 

This function is represented schematically by the dark strips in Fig. e and f, perpendicular to 
the direction of the projected wind vector Wxy. 



2.5. Residual Fluctuation Power 

The atmospheric fluctuation power remaining at the output of the instrument is determined 
by calculating the overlap integral between the unfiltered atmospheric power (Fig. ^) and the 
spatial and temporal filtering functions for the instrument as depicted in Fig. |3|. This is analogous 
to the application of a window function for assessing the response of an instrument to the CMB 
power spectrum, or to the use of the Optical Transfer Function for optical systems. Each of the 
three instrument configurations is considered below. 



2.5.1. Chopped beam 

The residual level of brightness temperature fiuctuations at the output of a chopped beam 
experiment after time averaging is given by 

/+00 r+oo 
/ {f^],y)&Q'^daxday (18) 
-OO <J — OO 

/-|~"00 p-\~oo 
^ I ^ {Ahli^isme)-^/Hc^l + al)-''/'}{s-^n\ne,^,^ax)} 

|sinc^(7r^av^^(ax^^^^: + ^y'Wy))^ daxduy. (19) 

The integral is dominated by the contribution from close to the origin, in which regime the 
spatial filter function can be approximated as Tr^^^hopO^.. The integral can then be reduced to the 
following expression: 

?;lt,rms^(ll-2cos2</,+ 16.7sin2</,)A/j2^(sine)-302j^„p^.;-^^^ (20) 
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where is the angle between the projected wind vector w^^y and the 9x axis. The expression 
should apply as long as the 3D turbulence model is a good approximation for (Tj^y) (Section p.2[ ). 
This requires that the chop angle ^chop ^ 2A/i sin e//iav and Wxyt^v ^ ^^av If either constraint 
is exceeded, equation ( pO[ ) represents an upper limit, since the regime of '2D turbulence' has 
reduced power at low angular wavenumbers. It is interesting to note that the rms level of the 
fluctuations goes down as t^^^^ , much more slowly than the usual tiJ^'^ . This is because the 
filtered fluctuations do not have a white noise spectrum; most of the power is concentrated at the 
low frequencies. 

It is clear that simple chopped observations are very susceptible to atmospheric fluctuations; 
in practice more complicated three and four beam chops have been used to alleviate this. 



2.5.2. Swept beam 



For a channel with Aacharu the atmospheric fluctuation power (T'J.y) can be considered 

approximately constant in the region of overlap between the spatial filter and the temporal filter 
(Fig. ^). For such a case, it can be shown that 



" out, rms 



{Tiy)G'Q'daxday 

-11/3- 
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(22) 



where a'^ is the angular wavenumber of the channel and (p is the angle that Wxy makes with 
the Ox-axis. The first factor is (T'gky) evaluated at the intersection of the channel spatial filter 
(shaded region in Fig. ^) and the temporal filter (dark strip in Fig. ^). The second factor is the 
taper imposed by the beam size of the aperture (a wide angle beam is insensitive to high angular 
wavenumbers). The third factor is the effective length of the temporal filter strip that overlaps 
with the channel. The fourth factor is the effective width of the temporal filter function. The 
third and fourth factors represent the area of overlap in Fig. |^; the first and second represent the 
level of fluctuations in the overlap region. The residual fluctuations are largest when the wind Wxy 
is blowing perpendicular to the sweep direction, i.e. (p = 90 degrees. 

The channels with a'x close to zero may have a significant response to the strong peak in 
(rj^y) near (a^,ay) = (0,0). Careful thought must be given to the tapering (or 'windowing') 
function applied over the swept slot on the sky, since this determines the profile of the channel 
as a function of ax. Little or no tapering leads to strong sidelobes on the channel response; too 
much tapering substantially reduces the resolution of the power spectrum. Specific cases can be 
computed numerically if the altitude /lav and thickness Ah of the layer are known. 
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2. 5. 3. Interferometer 

For an interferometer baseline B that is large compared to the diameter of the antennas, the 
circles representing the spatial filter response in Fig. ^ are widely spaced and the atmospheric 
power spectral density (Tj^y) can be considered constant over these regions. The overlap with 
the temporal filter (dark strip in Fig. ^) is maximum when the wind blows perpendicular to 
the baseline. For cases when the wind is close to perpendicular ((/> > 60°), the variance of the 
residual brightness temperature fluctuations at the time-averaged output of the instrument is 
approximated by 

p+oo p+oo ^ „ 



^out,rms = I I (^sky)G' Q da^^dtty (23) 

' B cos (/> \ ^ I 



-oo J — oo 



' D X -11/3 ~ 



exp 
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The first factor is (^sky) evaluated at the intersection of the channel spatial filter and the temporal 
filter (assumes '3D turbulence'). The second factor accounts for the taper of the primary beam. 
The third factor is the equivalent width of the Gaussian spatial filter along the direction of the 
temporal filter strip for both positive and negative ax- The fourth factor is the equivalent width 
of the temporal filter function. When the wind blows more parallel to the baseline there is very 
little overlap between the spatial and temporal filters (Fig. |3|f) and T^^^^_ j-ms very small compared 
to the perpendicular case. 

If the edge-to-edge separation AD of the apertures is small compared to their diameter d, the 
circles in Fig. ^ are close together, and it is necessary to make a numerical calculation to account 
for the change in in the region of overlap between the spatial and temporal filters. For the 

case where the wind is perpendicular to the baseline, equation (23) can be approximated by a one 
dimensional integral of the atmospheric power law and the spatial filter along the Ox direction, 
multiplied by the effective width of the temporal filter: 

rolt,rms ^ {AhliHsine)-'/'} |(^,,t,,)-i^| 5int, (25) 



where 



'int 



hOO 



a-^^/^G^dox. (26) 



This integral was evaluated for two different types of aperture: (1) a truncated Gaussian 
distribution of the electric field strength with the edge cut-off at the -10 dB level; (2) a truncated 
Bessel distribution of the electric field, cut off at the first zero. The first is typical for a dish 
illuminated by a feedhorn, and the second is obtained at the aperture of a corrugated horn. In 
each case, G was calculated by cross-correlating the E-field distribution of two apertures. It was 
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found that the result is almost independent of the transition from the 3D to the 2D turbulence 
regime. Also, for a given value of AD/d, the integral 5int oc (d/X)^^/^ (to a good approximation), 
i.e. S'int(d/A)^/^ is a function of AD/d. This is plotted in Fig. ^. At large values of AD/d (where 
AD ~ B), there is a —11/3 power-law dependence, as predicted by equation (p4[). 

2.6. The Unfiltered Fluctuation Strength 

The atmospheric brightness fluctuations at a given site and at a given time can be 
characterized by four parameters: the wind vector w, the altitude of the turbulent layer h^v, 
the fluctuation intensity A, and the thickness of the layer Ah. Of particular importance is the 
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combination AhJ, , which is the measure of fluctuation 'strength' relevant to the swept beam and 
interferometric observations of the CMB (eqs. [^] and |2^). 

In the next sections, data from the South Pole and Chile are analyzed to constrain these 
parameters. The models developed above are then used to estimate the sky brightness fluctuations 
that would be expected for an interferometer located at these sites. 

3. The South Pole 

The South Pole has been chosen as a site for several CMB anisotropy experiments over the 
past decade (see for example Meinhold & Lubin 1991, Tucker et al. 1993, Dragovan et al. 1994, 
Piatt et al. 1997). It is high (2800 m), extremely cold and dry, and is situated on an expansive 
ice sheet, with the surface wind dominated by weak katabatic airflow from higher terrain several 
hundred kilometers away to grid northeast (see discussion in King & Turner 1997). We characterize 
atmospheric brightness fluctuations at the South Pole using data from the Python telescope, a 
swept beam CMB experiment, obtained during the austral summer 1996-1997. 

3.1. The Python Experiment 

The Python telescope, in its configuration for the 1996-1997 season, employs a dual feed 
40 GHz HEMT based receiver. The receiver has two corrugated feeds separated by 2.75° on the 
sky, with separate RF chains, HEMT amplifiers, and backend signal processing. The post-detector 
output is AC coupled with a cut-off frequency of 1 Hz, and incorporates a low-pass anti-alias filter 
which attenuates the signal at frequencies above 100 Hz. The RF signal is separated into two 
frequency bands, centered at 39 GHz and 41.5 GHz with bandwidths of approximately 2 GHz and 
5 GHz, respectively. Data from the two bands are combined for atmospheric analysis. 

The feeds are in the focal plane of a 0.8 m off-axis paraboloidal mirror, resulting in two 1.1° 
beams on the sky which are swept through 10° at constant elevation at a rate of 5.1 Hz by a large 
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Fig. 4. — Plot of logio 'S'int(d/A)^/^ ^ ^ function of the ratio of the edge-to-edge separation to 
the diameter d of the aperture. The lower of the two lines is for a truncated Gaussian aperture 
distribution; the upper is for a corrugated horn. 



-16- 



vertical flat mirror. The two beams are at the same elevation and their sweeps partially overlap on 
the sky. Beam spillover is reflected to the sky by two sets of shields, one set fixed to the tracking 
telescope, and a set of larger stationary ground shields which also shield the telescope structure 
from the sun and any local sources of interference. Data are taken for ~ 30 s while sweeping 
and tracking a central position on the sky; the telescope is then slewed to another position a few 
degrees away. Between 5 and 13 pointings are stored in one data file, representing 5 to 10 minutes 
of observing time. Data were taken over 80% of the period from early December 1996 through 
early February 1997. 

3.2. Data Analysis Technique 

In order to differentiate atmospheric fluctuation power from instrument noise, the covariance 
of the data from the two beams is taken for the portion of each sweep in which their positions 
overlap on the sky, approximately 6°. Atmospheric brightness fluctuations are correlated between 
the two beams, while most of the instrument noise is uncorrelated. Thus the signal-to-noise of the 
correlated fluctuation power can be increased by averaging the covariance over many sweeps on the 
sky, allowing the atmospheric brightness fluctuation power to be estimated during stable periods 
when the system is receiver noise limited. The mean covariance represents the mean 'snapshot' 
fluctuation power in the 6° sweep, regardless of the number of sweeps that are subsequently 
averaged together. 

Several instrumental effects are accounted for when estimating the atmospheric brightness 
fluctuation power: 

1. The fluctuation power is corrected for the effect of the anti-alias filter roll-off and backend 
electronics delay. 

2. Correlated 60 Hz line noise is removed from the data. 

3. An offset dependent on the position of the sweeping mirror is correlated between the two 
beams. This is removed by subtracting the component of the signal that is constant on the 
sky over multiple pointings. 

We determine the effect of the 60 Hz line noise and stationary signal removal techniques on 
the true atmospheric signal by examining their effect on the data during periods when the data 
are dominated by atmospheric fluctuations. We calculate that the atmospheric power should 
be increased by 30% to compensate for the combined effect of these removal techniques and 
instrumental effects. There is an additional correlated signal due to the stationary ground shield 
and due to the CMB itself, which is not removed with these techniques. Therefore, the quartiles 
which we report should be taken as an upper limit of the true atmospheric signal. 
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3.3. South Pole Atmospheric Fluctuation Data 

/.From the Python data, we determine the South Pole brightness fluctuation power for a 
6° sweep at a mean elevation angle of 49° over 2 months of the austral summer (Fig. ^). The 
atmospheric fluctuations at the South Pole are bimodal in nature, with long periods of high 
stability broken by periods of high atmospheric fluctuation power. An examination of the 
meteorolgical records shows that the latter are always associated with at least partial cloud cover. 
The infrequent periods of high atmospheric fluctuation power are correlated with a (grid) westerly 
shift in the wind direction, and are weakly correlated with increased windspeed and precipitable 
water vapor content. These correlations indicate that the periods of high atmospheric fluctuation 
power are associated, at least in part, with synoptically forced moist air from West Antarctica, a 
condition which occurs infrequently at the pole (Hogan et al. 1982). Variability in the fluctuation 
power during periods of low atmospheric fluctuation power is consistent with instrument noise. 

From these data we construct a cumulative distribution function for the brightness fluctuation 
power (Fig. ^), and derive quartile values for brightness fluctuation power over the 2 month time 
period during which the data were taken. The cumulative distribution function for the data files, 
in which the covariance has been averaged for only a few minutes, is adequate for deriving the 50% 
and 75% quartile values. However, instrument noise dominates the distribution function at the 
25% quartile level; therefore data that are binned in 6 hour intervals are used to derive the 25% 
quartile value. The quartile values for the 2 months of the austral summer 1996-1997 are {0.20, 
0.51, 1.62} mK^, where the brackets are used to denote the 25%, 50% and 75% quartile values. To 
compensate for instrumental filtering effects described above, we increase these values by 30% to 
{0.27, 0.68, 2.15} mK^. The Python experiment was not operated during the austral winter, so 
data are not available for this period. However, precipitable water vapor and sky opacity quartile 
values are lower in the winter months (Chamberlin et al. 1997). It is therefore likely that the 
atmospheric stability improves during the austral winter as well. 

The sharp roll-off of the Python primary beam spatial filter (see Fig. ^) prevents an 
accurate determination of the underlying atmospheric angular power spectrum. The angular 
power spectrum for data that have not been tapered at the edges of the sweep exhibits power 
proportional to the inverse square of the angular wavenumber. This power law dependence is 
due to power from low angular wavenumbers (large angular scales) leaking into higher angular 
wavenumber channels. The observations demonstrate that it is desirable to taper CMB data to 
reduce contamination from the high atmospheric fluctuation power on large angular scales. 

3.4. Estimating the fluctuation intensity 

The value of Ah%^ (eq. [|l^) appropriate for the South Pole can be estimated from the rms of 
the Python measurements. Model atmospheric fluctuation power spectra for layers of turbulence 
seen at an elevation of 49° were computed using a full three dimensional integration (similar to 
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Fig. 5. — The brightness fluctuation power, in mK^, for 2 months of Python observations at the 
South Pole. The measured fluctuation power is for a 6° sweep on the sky, at a mean elevation angle 
of 49°. The data have been binned in 6 hour intervals to reduce instrument noise. In the lower plot, 
the scale of the vertical axis is enlarged about zero to show the fluctuation power during stable 
periods. Even after binning into 6 hour intervals, the variability in the fluctuation power during 
stable periods is consistent with instrument noise. 
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Fig. 6. — The cumulative distribution function of brightness fluctuation power for the Python 
data shown in Fig. |5[ Two distributions are shown, one for unbinned data files taken over a 
few minutes (solid line), and one for the data binned in 6 hour intervals (dash-dot line). The 
unbinned data files are adequate to determine the 50% and 75% quartile values, but the 6 hour 
binned cumulative distribution function is needed to achieve adequate signal-to-noise to estimate 
the 25% quartile value. The increased fraction of high fluctuation power data in the 6 hour 
binned cumulative distribution function results from binning brief periods of high fluctuation power 
together with periods of low fluctuation power. The quartile values are {0.20, 0.51, 1.62} mK^ 
(before compensating for instrumental effects). 
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Lay 1997). These were then multiphed by the Python spatial filter function (the sum of all but 
the DC channel) and integrated to give a value for AT^^^ j.^^. The model value for Ah^l^ was then 
scaled so that the model AT^^^^^^^ matched the measured value. Note that since we are estimating 
the level of rms fluctuation level from 'snaphots' of the atmosphere, there is no correction needed 
for time averaging, i.e. the temporal filter = 1 in equation ^. 

The result is a function of the ratio of the layer thickness to the average altitude of the 
layer, A/i//iav. For A/i//iav = 1, it was found that Ahli^ = 0.99AT^^^^^^; for Ah/h^^ = 2, 0.5 
and 0.25, the scale factor changed to 0.89, 1.13 and 1.36, respectively. Therefore the ratio does 
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not have a major effect on the inferred value of Ah^v ■ For Ah = h^v = 500 m, the measured 
quartile values give Ahti^ = 0.99AT^^^ = {0.27,0.67,2.13} mK^. Corresponding values for the 
quantity Ah^i^ / (uiK'^'ni) (which determines the residual level of fluctuations for swept beam and 
interferometer experiments - see eqs. and |2^) are tabulated in Table || for three different 
values of /lav These numbers are appropriate for an observing frequency of 40 GHz. The 
brightness temperature of fluctuations due to water vapor scales approximately as the square of 
the observing frequency (except close to the line centers at 22 GHz and 183 GHz), so that Ahav 
should be scaled as observing frequency to the fourth power. 





South Pole /lav 




Chile /lav 




Quartile 


500 m 1000 m 2000 m 


500 m 


1000 m 


2000 m 


25% 


150 340 800 


1500 


9400 


60,000 


50% 


370 840 2000 


7000 


44,000 


290,000 


75% 


1200 2700 6400 


29,000 


190,000 


1,200,000 



Table 1: Values of Ahli^ in units of mK^ m (see eqs. [^2| and p^) for the South Pole as estimated 
from Python V data (section ^) , and for the Atacama Desert in Chile estimated from phase monitor 
data (section ^). The altitude /lav of the turbulent layer is not known. The Chile numbers have 
a 50% uncertainty associated with the conversion from refractive index to brightness temperature. 
Values are appropriate for 40 GHz, but can be scaled to other frequencies based on the emissivity 
spectrum of water vapor. 



3.5. Altitude of the fluctuations 

By combining the Python data with radiosonde wind measurements, it is possible to determine 
the altitude of the fluctuations during periods of bad weather. This is illustrated by Fig. 0, 
which shows the measured emission as a function of angular position 6x on the sky and time 
t. The plot represents an interval of 30 s, at a time when the wind was blowing parallel to the 
sweep direction. The stripes are produced by blobs of water vapor moving from left to right; the 
diagram shows that a blob moves through an angle of 7° in about 13 s. A 7° angular distance 
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corresponds to a physical length of 0.12/iav/ sin(49°), where /lav is the average altitude of the 
fluctuation. The radiosonde launched 2 hours after this dataset indicated a fairly uniform wind 
speed of w =16 m s~^ for the lower 2 km, so a blob is expected to move 208 m in 13 s. Solving for 
/lav gives 1300 m. The slope of the stripes is proportional to hi^^/w. 

In order to average the data together from many 30-s intervals, we compute the power 
spectrum for each time-angle plot and average those together. The power spectrum for a 30-s 
interval is calculated by computing the Fast Fourier Transform of the time-angle data from 
each of the two feedhorns (using a Hann taper to minimize sidelobes), and then calculating the 
covariance between the two transformed datasets. Figure 0b shows the power spectrum averaged 
over one hour of data (containing the 30-s interval shown in Fig. ^). The power is distributed 
along a radial line, perpendicular to the striping. The gradient of this line is proportional to 
if/Ziav; the overlaid radial lines represent (from vertical) altitudes of 0, 500, 1000, 1500, 2000, 2500 
and 3000 m, calibrated using w =16 m s~^. The fluctuations have h^v ~ 1300 m. The physical 
periodicity is given by w/i'; e.g. v = 0.1 Hz corresponds to a fluctuation with period 160 m. 
There is little power present at the origin because the DC level was removed from the time-angle 
plane before the Hann taper and Fourier transform were applied. The contours appear to fall 
off faster than would be expected for Kolmogorov turbulence. This is partly due to the effect of 
the primary beam taper in the angular wavenumber direction, which can be represented by a 
Gaussian centered on zero with a FWHM size of about 33 rad~^, but may also indicate that these 
bad weather fluctuations do not follow a Kolmogorov power law. 

Average power spectra were computed for all the hours when the wind was parallel to the 
sweep direction (approximately every 12 hours). Two more examples during periods of bad 
weather are shown in Fig. ^ and d. The first indicates fluctuations at an altitude of ~ 500 m; the 
second shows two components: one at ~ 300 m and another at much higher altitude (> 3 km). 
Unfortunately it was not possible to detect structure in the power spectrum during stable periods; 
the emission is too weak. 

In most cases that were measured, the altitude determined for the strong fluctuations, which 
varies from 300 m to well over 3 km, agreed well with the altitude at which the relative humidity 
was a maximum (measured by radiosonde launches). This strengthens the case for the connection 
between clouds and strong fluctuations, and may explain the possible non-Kolmogorov nature of 
the power spectrum. In the cases of high altitude turbulence (/lav > 3 km), however, there was no 
corresponding maximum in the relative humidity and there is generally little water vapor present 
in the atmosphere. Another mechanism must be at work in these cases. 

4. The Atacama Desert in Chile 

The Atacama Desert in Northern Chile is extremely dry, and is the proposed location for the 
next generation of large millimeter-wave arrays, as well as the Cosmic Background Interferometer 
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Fig. 7. — a) Time-angle plot for 30 s of Python V data during bad weather, illustrating the striping 
that results when the wind blows blobs of water vapor along the sweep direction, b) Power spectrum 
for a 1 hour period which includes the sample shown in (a). Contours are relative to maximum 
(0.9, 0.8, 0.7, 0.2, 0.1, 0.09, 0.08, etc.). Note that the power is concentrated in a line that 
is perpendicular to the striping in (a). The radial lines correspond to different altitudes, based 
on the windspeed measured by radiosonde launches. The y-axis corresponds to ground level, and 
subsequently shallower lines represent 500, 1000, 1500, 2000, 2500 and 3000 m, respectively. The 
fluctuations in (a) are at roughly 1500 m above ground, c) Another example showing fluctuations 
at less than 500 m. d) Two fluctuation components; one very low, the other very high. 
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experiment. The Mobile Anisotropy Telescope has been deployed there for two seasons during 
1997-1998 (Torbet et al. 1999, Miller et al. 1999). Monitoring of the atmospheric stability at Cerro 
Chajnantor (the proposed site for the Atacama Large Millimeter Array) has been underway for 
over 4 years, using a site test interferometer (Radford et al. 1996). 

4.1. The site test interferometer 

The interferometer used in Chile consists of two dishes 1.8 m in diameter, separated by 
an East- West baseline of 300 m, that observe the 11 GHz CW beacon from a geostationary 
communications satellite at an elevation of 36° and an azimuth of 65°. The instrument measures 
the phase difference between the two received signals, which depends on the difference in the 
electrical path lengths along the two lines of sight to the satellite. The random component of the 
fluctuations is dominated by the non-uniform distribution of water vapor in the troposphere being 
blown over the interferometer, since water vapor has a high refractive index compared to dry air. 
Quartile values for the rms difference in path length over 1996 are {136, 300, 624} /xm. The data 
processing is described by Holdaway et al. (1995). These values have not been scaled to zenith. 

Water vapor is the principal source of both path length fluctuations and brightness 
temperature fluctuations. We wish to estimate the latter from a measurement of the former, 
which requires a conversion factor. This is described in § [4.2| . It is also necessary to determine 
A (the intensity of fluctuations integrated through the atmosphere) from the interferometer 
measurements. This requires a model, and is described in § [4.3| . 

4.2. Converting rms path to brightness temperature 

The millimeter and submillimeter absorption spectrum of water vapor is dominated by a 
series of rotational line transitions, the lowest of which are at 22 GHz and 183 GHz. At frequencies 
close to the line centers, the optical depth r for a given amount of precipitable water vapor 
(PWV: the depth of liquid water obtained if all the vapor is condensed) is well-determined, but 
between the lines experiments have shown that the absorption is higher than expected from theory 
(e.g. Waters 1976, Sutton and Hueckstaedt 1996). The reason for this is still not clear, although 
several hypotheses have been advanced. At 40 GHz the theoretical absorption from lines is 0.012 
per centimeter PWV. The Gaut and Reifenstein (1971) model for the 'continuum' absorption 
increases this to 0.03 cm~^, and is believed to be a reasonable estimate at the frequencies of 
interest here (Sutton and Hueckstaedt 1996). If the physical temperature of the atmosphere is 
250 K, then 1 cm PWV has a brightness temperature of 0.03 x 250 K= 7.5 K. The excess path 
resulting from 1 cm of water vapor at 250 K is approximately 7 cm (Thompson, Moran & Swenson 
1994), so that 1 cm of excess path corresponds to a brightness temperature of 1.1 K, and 1 mK 
brightness corresponds to 9 /um excess path. This conversion should be considered an estimate 
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with an uncertainty of order 50%. A hard upper hmit of 22 /xm mK ^ is set by the contribution 
from hne emission only (no continuum term). 

Using the value of 9 /um mK~^, the rms path fluctuation quartile values are mapped to 
brightness temperature fluctuations of {16, 34, 69} mK at 40 GHz. 



4.3. Estimating A 

Since the phase monitor measures the path difference between two lines of sight, its spatial 
filtering properties are analogous to the chopping instrument depicted in Fig. |^ and d. The main 
difference is that the lines of sight are parallel to one another through the atmosphere, rather 



than diverging from the point of observation. As discussed in §2.5, calculation of the residual 
fluctuation power after spatial and temporal filtering for this configuration requires a model that 
includes the thickness of the layer containing the fluctuations and the detailed geometry of the 
observations. The relevant analysis can be found in Lay (1997). 

The input model parameters are: baseline 300 m East- West; elevation 36°; azimuth 65°; 
layer thickness 500 m. The windspeed and altitude of the layer are not needed for calculation of 
A. To obtain rms brightness temperatures of {16,34,69} mK with the above model parameters 
requires that A = {9.4 x 10~^,4.4 x 10~'^,1.9 x 10~^} mK^ m"^/^. The layer thickness 
has only a small effect on the calculation; if instead the layer is actually 2 km thick, then 
A = {8.0 X 10-^ 3.7 X 10^^, 1.6 x lO^^} mK^ m~^/^. 

The average altitude of the turbulent layer is much more important, but is not known. The 
quantity Ahti^ / (mK'^m) is tabulated in Table [l| for different values of /lav and Ah = 500 m. This 



is the value that determines the residual level of fluctuations in a CMB experiment (eq. |24|) 



5. Example 

The analysis and results presented above are applied to calculate the contribution from 
atmospheric brightness fluctuations to the noise measured by a small interferometer, of the kind 
that might be used for measurements of the CMB anisotropy. 

We specify a field of view with FWHM of 6'b = 3° = 0.053 rad at a frequency of 40 GHz 
(A = 0.75 cm). The diameter of the aperture required for a Gaussian distribution truncated 
at the —10 dB level is d = 1.15A/0b = 16.3 cm. For a corrugated horn, the diameter is 
d = 1.32A/0b = 18.7 cm, and it is this case that we now investigate further. 

The total thermal noise (for the combined sine and cosine components, no atmospheric 
fluctuations) for a 2-element interferometer with system temperature Tgyg and bandwidth Az^ is 
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given by (see Thompson, Moran & Swenson 1994) 

T.,™,™ = 0.28 (1^) {^y ' (t)-"^ OK. (27) 

The rms atmospheric contribution is given by the square root of equation ([25|): 

Tout,r.s={Ahl(^f\sine)-'y'w-y't-y^ (28) 
The value of S'int is determined from Figure H, and depends on the separation of the apertures. The 
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value of Ahav depends on the conditions at the site, as shown in Table as does the projected 
windspeed Wxy 

An ideal site should not significantly compromise the sensitivity of the experiment. We can 
determine the range of Ahti'^ for which Tout,rms < ^therm,rms- From Fig. ^, the maximum value 
of \og Sint{d/ X)^^^ is approximately +0.3, which apphes when the wind is perpendicular to the 
baseline for an interferometer with corrugated horns that are very close together. In this example, 
d/X = 24.9, which implies S'int < 3.8 x 10~^. Substituting these values into equation (|28|), along 
with a representative elevation of 60° and a projected windspeed of 10 m s~^, we obtain 

Tout^ms < 8.0 X 10-3 (a/.«/3)'/%-,V2^i/2^-i/2_ (29) 

For this to be less than the thermal noise contribution we require Ahli^ < 1.2 X 10^ mk^ m . 
Reference to Table |l] indicates that the South Pole should satisfy this requirement most of the 
time during the summer months, unless the fluctuations are higher than 2 km above the ground. 
At the Chile site, it appears that the instrument noise for this configuration would be dominated 
by atmospheric fluctuation power, regardless of the altitude of the fluctuations. It is important, 
however, to realize that this example was calculated for the extreme case of the wind blowing 
perpendicular to a baseline for which the apertures are almost touching. The situation is improved 
by separating the apertures to observe smaller angular scales (Fig. ^), or when the wind blows 
parallel to the baseline (Fig. ^). 



6. Summary 

1. The impact of non-uniform emission from the atmosphere on measurements of the CMB 
is assessed. Chopped beam, swept beam and interferometric measurement schemes are 
analyzed. 

2. The analysis is based on a model where the atmospheric fluctuations are confined to a 
turbulent layer. Data from phase monitors indicate that there is a transition from a three- 
to a two-dimensional regime on scales comparable with the thickness of the layer. The 
outer correlation scale of fluctuations is much larger than the layer thickness. Previous 
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analyses have assumed a much smaller outer scale length and therefore underestimate the 
level of fluctuations on large scales. The impact of atmospheric fluctuations is assessed by 
considering the instruments as a combination of spatial and temporal filters that act on the 
power spectrum of the turbulence. 

3. Data from the Python V experiment during summer at the South Pole are analyzed to 
determine the level of fluctuations. The distribution is bimodal, with stable conditions for 
75% of the time and much stronger fluctuations 25% of the time. 

4. During normal stable conditions the rms brightness temperature variation across a 6° strip 
was less than 1 mK, at a frequency of 40 GHz. It was not possible to determine either the 
average altitude or the power spectrum of these very weak fluctuations. 

5. The bad weather fluctuations have rms brightness temperature variations of up to 100 mK 
across a 6° strip of sky, at a frequency of 40 GHz. They are only present when there is at 
least partial cloud cover. The windspeed as a function of altitude was used to determine the 
average altitude of the fluctuations. This varied from 500 m to over 4 km, and in most cases 
corresponded to the altitude at which the relative humidity was a maximum. The power 
spectrum measured for these bad weather fluctuations appears to fall off faster than the 
predicted Kolmogorov power law, although this may be due to the taper imposed by the 
primary beam. We conclude that these strong fluctuations are probably associated in some 
way with cloud activity. 

6. Path fluctuation data from a satellite phase monitor located in the Atacama Desert in Chile 
were used to estimate the corresponding level of brightness temperature fluctuations at 

40 GHz for that site. Although this conversion is uncertain by a factor of 2, and the result 
depends strongly on the unknown altitude of the turbulent layer, the two months of South 
Pole data indicate a significantly lower level of fluctuations compared to the Chile site, as 
shown in Table |^. 

7. The theoretical analysis and experimental measurements were combined to predict the 
residual atmospheric noise that would be present at the output of a small interferometer 
with antennas of diameter ~ 25A that could be used to measure CMB anisotropy. It was 
found that the atmospheric contribution is likely to be small compared to the thermal noise 
when the apertures are well-separated, but that the atmosphere can dominate in cases where 
the edges of the apertures are separated by only a few wavelengths. For the latter case, a 
good site, such as the South Pole is critical. 

8. The analysis shows that the average altitude of the turbulent layer has a big impact on 
the suitability of a site for CMB measurements. This parameter is not well constrained by 
the data from either site. Measurement of this altitude should be a priority for future site 
testing experiments. 
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